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Abstract 

Nonlinear optical properties of extracted higher fullerenes - C70, C76, C78, and Cs4 - are theoretically 



P^h ' investigated. Magnitudes of off-resonant third-harmonic-generation are calculated by the intermediate exciton 

£SJ ■ theory. We find that optical nonlinearities of higher fullerenes arc a few times larger than those of Ceo- The 

magnitudes of nonlinearity tend to increase as the optical gap decreases in higher fullerenes. 
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It has been found that Ceo thin films show large optical nonlinearities. 1-4 ' The observation is attractive 
in view of scientific interests as well as technological applications. The magnitude of the third-order nonlinear 
susceptibility, Xthg( w ) = X^(3w;w,w,w), for third harmonic generation (THG) is of the order of 10~ 12 esu 
to 1CP 11 csu. This large response is comparable to the values measured in polydiacetylenes. The optical 
spectra of C70 4 ' and higher fullerenes (C76, C78, Cg4, etc.) 5 ' 6 - 1 have also been obtained. In order to explain the 
interesting experiments, several theoretical investigations 7-11 ' have been performed. We have applied a tight 
binding model 7 -* to Ceo, and have analyzed the nonlinear optical properties. Coulomb interaction effects on 
the absorption spectra and the optical nonlinearity have been also studied. 10 ' We have found that the linear 
absorption spectra of Ceo and C70 are well explained by the Frenkel exciton picture 11 ' except for the charge 
transfer exciton feature around the excitation energy 2.8 eV of the Ceo solids. 12 ' Coulomb interaction effects 
reduce the magnitude of the optical nonlinearity from that of the free electron calculation, 10 ' and we have 
discussed that the local field enhancement might be effective in solids. 

In the previous paper, 13 ' we have investigated geometric effects on optical properties in higher fullerenes. 
We have obtained the optical absorption spectra at a certain combination of the pentagonal carbons, by using 
projected wavefunctions onto selected pentagons in order to calculate dipole moments. The contributions from 
a part of the fullcrcnc to the optical spectra have been extracted. We have found that the optical excitations 
in the energy region smaller than about 4 eV have most of their amplitudes at the pentagonal carbons. The 
oscillator strengths of absorption projected onto these carbons almost accord with those of the total absorption. 
We have also found that the spectral shapes of the total absorption are mainly determined by the geometrical 
distributions of the pentagons in the fullcrcnc structures. 

The main purpose of this paper is to investigate nonlinear optical properties in higher fullerenes, further. 
The calculation method of the THG has been used for Ceo in the previous work. 10 ' Recently, a calculation of 
the THG in isomers of C78 by a free electron model has been reported. 14 ' However, the Coulomb interaction 
effects, whose importance we have found in Cgo, have not been discussed for higher fullerenes. The present 
paper gives rise to a new contribution to this subject. We focus on the off-resonant THG in order to estimate 
possible magnitudes of the nonlincarities of each isomer. The Coulomb interaction strengths are also changed 
in a reasonable range, because realistic strengths are not well known in higher fullerenes. 

In the present paper, the carbon network of the fullerene surface is taken into account by the hopping integral 
t between nearest neighbor sites. The Coulomb interactions are taken into account by the parametrized Ohno 



potential, W(r) = 1/ yJ(l/U) 2 + (r/roV) 2 , between two electrons with distance r. Here, U is the interaction 
strength at the distance r — 0, V is the strength of the long range part, and r$ the mean bond length. The 
Coulomb interaction is treated by the restricted Hartree-Fock approximation and the intermediate treatment of 
excitons. 11 - 1 The THG is calculated by the sum-over-states method which has been used in ref. 7. In calculating 
the expectectation values of the dipolc moment, the lattice coordinates contained in the geometry package of 
higher fullerenes 15 -* are used. We will discuss properties of the THG in seven kinds of the extracted fullcrcnc 
isomers: C70, C76, C78, and Cs4- From our experiences of investigation on optical properties, 10 ' 11 ' we can 
assume V — U/2 in the following. The onsitc Coulomb strength is varied within < U < At. 

Figure 1 shows the absolute value of the off-resonant THG |Xthg(^)I plotted against the Coulomb interaction 
strength U. The different plots indicate diffcrenct kinds of isomers. The four isomers - C70, -D2-C76, D3-C78, and 
one kind of C2V-C78 [C^v by Kikuchi's notation (ref. 16)] - exhibit similar magnitudes of optical nonlincarities. 
On the other hand, the other three isomers - one kind of C2 V -Crs [C-2v by Kikuchi's notation (ref. 16)], £ ) 2d-Cs4, 
and -D2-C84 - show larger optical nonlinearities than those of the former isomers. This is mainly due to the 
smaller energy gap of the latter isomers, even though the negative correlation between the THG and the energy 
gap is not so complete through all the isomers. The decrease of THG from the free electron model (U = 0) 
to the case with U = At is by the factor about 0.1 for all the isomers. The similar fact has been found in 
the calculation of Ceo- This would be a general property in various kinds of higher fullercncs. The overall 
magnitudes of the THG are around 10 -12 esu for most of the Coulomb interactions considered. 

In Fig. 2, the relations between the absolute value of the off-resonant THG IXthgWI ano - t ne energy gap 
are shown for three Coulomb interaction strengths: U = Oi, 2i, and At. Here, the energy gap is defined as the 
optical excitation energy of the lowest dipolc allowed state. This is called the optical gap, in other words. For 
each Coulomb interaction, the seven plots cluster in a bunch. When the energy gap becomes larger, the THG 
tends to decrease. However, the correlation between the THG and the energy gap is far from that of a smooth 
function. The correlation is merely a kind of tendency. Therefore, the decrease of the energy gap in higher 
fullercncs is one of origins of the larger optical nonlinearities of the systems. The actual magnitudes would also 
be influenced by detailed electronic structures of isomers. 

In the calculations of Cgo, the magnitudes of the THG at the energy zero are about 1 x 10 12 esu in the free 
electron model (ref. 7), and about 1 x 10~ 13 esu for U = At and V — It (ref. 10). In the present calculations 
of higher fullercncs, the magnitudes are a few times larger than those of C@o- Thus, this paper predicts that 



nonlinear optical responses in higher fullerenes are generally larger than in Ceo- In the previous paper, 10 ) we 
have discussed that the local field correction factor is of the order 10 in Ceo solids. As the distance between the 
surfaces of neighboring fullerene molecules in C70 and C76 solids is nearly the same as in Ceo solids, we expect 
that local field enhancement in thin films of higher fullerenes is of the similar magnitudes as in Cgo systems. 

In summary, we have investigated nonlinear optical properties of higher fullerenes. Theoretical off- resonant 
THG has been calculated by using the intermediate exciton theory. We have found the optical nonlinearities 
of higher fullerenes which are larger than in Ceo- The magnitudes of THG tend to increase as the optical gap 
decreases in higher fullerenes. 
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Figure Captions 



Fig. 1. The absolute value of the off-resonant THG |Xthg(0) plotted against the Coulomb interaction strength 
U. The closed and open squares are for C70 and D2-C76. The closed circles are for D3-C78, and the open 
and crossed circles are for two kinds of C2«-C78. The closed and open triangles are for D2d-Cs4 and D2-C84, 
respectively. 

{3} 
Fig. 2. The absolute value of the off-resonant THG |Xthg(0)I m seven isomers of higher fullercncs, plotted 

against the energy gap (shown in units of t). The squares, circles, and triangles are for U = Ot, 2t, and At, 

respectively. The left axis is in the logarithmic scale. 



